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Abstract 

We consider those Gaussian Unitary Ensembles where the eigenvalues have pre- 
scribed multiplicities, and obtain joint probability density for the eigenvalues. In the 
simplest case where there is only one multiple eigenvalue t , this leads to orthogonal 
polynomials with the Hermite weight perturbed by a factor that has a multiple zero at 
t. We show through a pair of ladder operators, that the diagonal recurrence coefficients 
satisfy a particular Painleve IV equation for any real multiplicity. If the multiplicity 
is even they are expressed in terms of the generalized Hermite polynomials, with t as 
the independent variable. 
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1 Introduction 



Random matrix ensembles originally conceived to explain the statistical properties of the 
energy levels in heavy nuclei jHOl has recently seen applications in transport in disordered 
systems, string theory and various areas of pure and applied mathematics. In addition to 
classical quantities of interest such as the correlation functions, the average of the product 
of characteristic polynomials of random matrices were under investigation starting from 
Brezin-Hikami paper |3] (see also j3] and references therein). 

From the Painleve equations point of view the average of a power of characteristic poly- 
nomial in Gaussian unitary ensemble gives a r -function of the rational solution of Painleve 
IV equation. For the integer powers this can be seen from the original Kajiwara-Ohta de- 
terminant formula for the rational solutions of PIV ( JHl; c.f. [22]) and it was later explored 
by Forrester and Witte [T3] . 

In this paper we consider the degenerate gaussian unitary ensembles. That is we restrict 
ourselves to the nonlinear subspace of Hermitian matrices having prescribed spectrum de- 
generacy. Various statistical properties on eigenvalues of such matrices can be asked. The 
first natural question we answer is the determination of the joint probability density of the 
eigenvalues when these have some multiplicity. It happens that as in the case of classical en- 
sembles (see 20] ) the joint probability density has the form of product of pairwise differences 
between the different eigenvalues taken in the powers depending on the multiplicities. Thus 
we naturally arrive to considering orthogonal polynomials with Hermite weight perturbed 
(multiplied) by a product of linear factors. 

These types of weights also appear in the random matrix theory in consideration of 
the averages of characteristic polynomials (see j2]), although in that case the zeroes 
of these factors are the external variables to the matrices of the ensembles. We note that 
an orthogonal circular random matrix ensemble with fixed degenerate eigenvalue at 1 was 
considered by Snaith in j^E] in conjectural relation to number theoretic questions on L- 
functions of elliptic curves. More general Jacobi circular ensembles were studied recently in 
jTTj. A general approach to the joint probability density of ensembles of various type was 
suggested recently in p. 

We are also motivated by the theory of Calogero-Moser-Sutherland systems. The ground 
states of these systems at appropriate interaction parameter coincide with the joint prob- 
ability densities of eigenvalues in the classical ensembles. The joint probability density 
for degenerate ensembles coincides with the factorized wave function for the appropriate 
multi-species generalisation of Calogero-Moser problem considered in This type of 

generalisation is integrable in the case of two types of particles p], [21] ■ More remarkably, 
Sergeev and Veselov showed that the corresponding quantum Hamiltonian can be obtained 
by applying a restriction procedure on the Calogero-Moser-Sutherland Hamiltonian in the 
infinite dimensional space to the appropriate discriminant [22] ■ We plan to elaborate these 
relations in future. 

In the context of orthogonal polynomials, perturbations of the standard weights such as 
the Jacobi weight by special factors is an important topic of investigation, where the problem 
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is the determination of the recurrence coefficients from the weights (see [21] and the 
references therein). In particular, it was noted by Magnus in |18J that often the variations 
lead to the recurrence coefficients which are solutions to the nonlinear equations. In some 
cases the appearance of Painleve IV for the certain exponential weights was established ^H] • 
More recently, it was shown in [TU| that the diagonal recurrence coefficient associated with 
the Hermite weight perturbed by special discontinuous factor satisfies a particular Painleve 
IV. 

In this paper we show that when the Hermite weight is perturbed by a linear factor having 
multiple zero the diagonal recurrence coefficients satisfy a particular two parameter Painleve 
IV equation. This property in fact holds for an arbitrary real power of the linear factor. Our 
approach is direct, it is based on an extension of the ladder operators technique developed 
in In Section 3 we describe this method, suitable for orthogonal polynomials where 
the weight has isolated zeros, in particular we derive a pair of fundamental compatibility 
conditions (Si) and (>S'2) . In Section 4 we make use of these to generate non-linear difference 
equations satisfied by the recurrence coefficients. These difference equations when combined 
with the Toda equations give a PIV equation satisfied by the recurrence coefficients • 

In the cases of the weights arising from degenerate gaussian ensembles having one multiple 
eigenvalue t of degeneracy K the multiplicity of the linear factor in the weight is 2K . In 
this case the recurrence coefficients are rational as functions of t . The theory of rational 
solutions to PIV [221 results in the expression of the recurrence coefficients through the 
generalized Hermite polynomials. 

We also mention here that the Hankel determinants associated to the Hermite weights 
perturbed by a factor are related to the Hankel determinants of the Hermite weight with the 
addition of 5 -function and its derivatives For an alternative derivation using Heine's 
multiple integral see [S]. 



2 Non-generic random matrices 

Let T-Ln be the space of Hermitian matrices of size N and let m = {mi, 1712, . . . yfrLk) be a 
partition of N . Consider the (nonlinear) subspace Ti^ in consisting of matrices having 
the eigenvalues with prescribed multiplicities mi, ... , . That is we suppose the spectrum 
{Ai, . . . , Aat} of an arbitrary element A G H"^ has the multiplicities described below: 

/il = Ai = . . . = Ami, 

= Ami+1 = . . . = Xmi+m2i 



fJ^k — Ami+m2+...+mfe_i+l — ... — Aat, (2-1) 

where we have renamed the eigenvalues as /ii, . . . without repetitions. 
As every Hermitian matrix A is diagonalizable, we have 

A = UAU-\ (2.2) 
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where A = diag(Ai, . . . , Atv) , and U is unitary. The matrix U is constructed out of a 
certain orthonormal basis where A becomes diagonal. Such a basis is defined up to unitary 
transformations leaving the eigenspaces invariant. Therefore U is determined as an element 
of the homogenous space 



where the direct product U{mi) x . . . x [/{rrik) of unitary matrices of orders mi, . . . , is 
embedded into U{N) as diagonal block. More precisely, in order to determine U uniquely, 
we also assume that the eigenvalues /ij, /ij having equal multiplicities rrii = mj are such 
that Hi < fij if i < j . 

Although the subspace H^, of Hn is a measure zero set, we may nonetheless construct 
a natural probability measure of the matrices lying in it. The metric 



is well-defined in the subspace Ti.^ . Therefore this metric also naturally defines a mea- 
sure on the subspace Ti"^ , via the Riemann volume formula. It happens, just like in the 
case of Hermitian matrices with distinct eigenvalues, with the spectral decomposition ()2.2j) . 
the measure on Ti.^ is a product of a measure on the eigenvalues and a measure on the 
homogeneous space ()2.3|) . 

Proposition 1 The metric Ii2.4\ ) restricted to the subspace H]} has the form 



U e U{N)/U{mi) X ... X f/(mfc). 





(2.4) 



k 



{dsf 



^ niid^il + 2 ^ (/ii - Hjf{dsijf, 



where 



{ds,,y = {U-'dU)^p{U 'dU)^fs. 




(2.5) 



The corresponding volume form on is 



k 




(2.6) 



where dv{U) is invariant measure on the homogeneous space \2.^) . 



Proof. From decomposition ()2.2|) we obtain 



dA = U {dA + U*dUA - AU*dU) U*. 



Then the metric (j2.4j) can be rewritten as follows: 



{dsf = tr {{dAf + 2{5UA - ASU)dA + {A5Uf + {5UAf - 25UA^SU) , 
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where 6U := U ^dU and we have used the cychc property of the trace. Simphfying this 
further we arrive at 

{dsf = tr j {dAf + 2 J] {XiXjSUijSUji - X'^i5Uij5Uji) 
and since 6U is anti-Hermitian we get (c.f., e.g., [Hj, 123) that the above reduces to 

N 

i=l l^i<j^N 

RecaUing the degeneracy conditions (2.1) we note that some of the terms vanish and the 
restricted metric takes the form 

k 

(c/s)2 = ^m,(rf/ii)2 + 2 J2 {^^^- ^^j?{dsij)\ (2.7) 

i=l l^i<j^k 

where dsij is defined in (|2.5p . The second sum in (|2.7p is well defined in the homogeneous 
space. 

To determine the corresponding measure we fix locally the section of the representatives 
of the coset classes and consider coordinates Ua/s such that dua/s = (U^^dU)ai3 , where 
the indices a < (3 are such that {a(3) ^ A. Here A is the diagonal block containing 
U{mi) X ... X U{mk) . Such local coordinates u^p will exist if the section is chosen to 
satisfy 5Uij = when (ij) G A . Then taking the real and imaginary parts ^Uap, as 
real coordinates the metric (j2.7|) becomes a diagonal metric gu and the term 2(/ij — /Xj)^ 
appears 2mimj times along the diagonal. From the Riemann volume formula, the measure 
corresponding to ()2.7|) is 

M k 
i=l i=l (q0)?A 

a<l3 

where 

M k 

1/2 



nv^= n 2"»-^(/i,-/i,)''"'"^^n^ 

1=1 ^<i<j<k i=l 

and M = dimTC^ . Thus we obtain the result (|2.6|) with the measure dh'{U) given by 

k 

du{U) = l[my' H {U-'dU)o^p(U~'dU)o.p. 

i=l (Q^f)^A 
a<l3 

Remark 1. A large class of generalized random matrix ensembles was recently consid- 
ered in [l] where a formula for joint probability density of the eigenvalues was obtained. 
Expression ()2.fi|l may be obtained from that work. 
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Remark 2. One way to generalize Proposition 1 is to consider real symmetric matrices 
with multiple spectrum, then angular variables are given by a factor in the orthogonal group. 
Same arguments as above lead to the following joint probabihty density of eigenvalues 

k 

l^i<j^k i=l 

Another possibility is to consider degenerate circular ensembles, that is unitary (or other) 
ensembles with given spectrum multiplicities. In this case the calculation of joint probability 
density results in taking the appropriate powers of nontrivial Cartan roots. 

It is a well-known result of Random Matrix theory that the partition function of any 
unitary invariant matrix ensemble defined by the multiple integral, 

Ar,[w] := — ... Yl {xi-XjfY[w{xk)dxk, (2.8) 

■ l<i<j<N k=l 

has the alternative representations, namely 

AjvM = det( f x'-^^w{x)dx^ (2.9) 

= det I / Pi{x)pj{x)w{x)dx j , (2-10) 

\J a J i,j=0 

where pi{x) is an arbitrary monic polynomials of exact degree /. Now if we orthogonalise 
these with respect to the weight w over [a,b], namely, 

b 

Pi{x)pj{x)w{x)dx = hiSij, 



a 



where hi, i eN is the square of the norm, then ()2.8|) becomes, 

N-l 

AtvM = n (2-11) 

j=0 

For the generic Gaussian Unitary Ensembles, w{x) = exp(— x^), x G M. In the case of 
a single degenerate eigenvalue t with K fold degeneracy and the rest, n eigenvalues are 
distinct, such that N = n + K , we find, by relabeling, /ii = t, fi2 = Xi, . . . , fik = Xn, the 
partition function reads. 



/oo 
e-^'' D^{t)dt, (2.12) 
-oo 

^"^^^ = d / / n - ^^•)' n(^' - ty^'e-^'dxi. (2.13) 



where 



l<i<j<n 
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We note the partition function expressions ()2.12|) . ()2.8p are defined here up to constant 
muhiples that come from the integration over the corresponding homogeneous spaces. 

The weight of orthogonal polynomials associated with integral (|2.13|) is the Hermite 
weight multiplied by an isolated zero, that is, 

t) = exp(— — tp"^, x, tGM. 

Other crucial characteristics of Random Matrix ensembles are the correlation functions 
of the eigenvalues. These are obtained by calculating the partition function type integrals 
()2.8p when some of the eigenvalues are fixed. In the case of single degenerate eigenvalue 
those correlation functions that involve the multiple eigenvalue coincide with the averages 
of the powers of characteristic polynomial for the appropriate standard Gaussian unitary 
ensemble, as it is immediately seen from (j2.12p - (j2.13|) . These averages were obtained in the 
determinant form in [5]. 

3 Ladder operators 

We now develop a differentiation formula for the polynomials Pn{x) orthogonal with respect 
to the weight wo{x)\x — t\'^ on the real line, for any smooth reference weight wq and for 
general 7 > . The derivation given here is similar to what was previously known [3 IH] , 
but adapted to the situation where the weight vanishes at one point. 

From the orthogonality condition, there follows the recurrence relations; 

ZPn{z) = Pn+l{z) + anPn{z) + l3nPn-l{z), 

with the initial conditions Pq{z) = 1, and Pop_i{z) = 0. The diagonal recurrence coefficients 
can then be expressed as 

a„ = pi(n) - pi(n + 1) (3.1) 
where Pi(^) are defined by expansions 

p^{z)=z^ + p,{n)z^-' + ... (3.2) 

The coefficients of the orthogonal polynomials will also have t dependence due to the t 
dependence of the weight although we denote the polynomials as Pn{z) . 

Since Pn{z) is a polynomial of degree n, its derivative is a polynomial of degree n — 1 
and can therefore be expressed as a linear combination of pk{z), k = 0,1, ...,n — 1 , namely, 

n-l 

P'niz) = ^Cn,kPk{z). (3.3) 
fc=0 

To determine the coefficients Cn,k we use orthogonality relations and the formula 

I "/"It 

d,\x - = s{x - t){{x - ty - (t - xy) + 7 '^" ' . (3.4) 
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We have 

h 



Cn,k = / Pniy)Pk{y)wo{y)\y -t\^dy 



oo 

oo 



/ Pn{y)Pk{y){wQ{y)\y -t\^ + Wo{y)dy\y -t\^)dy 



hi 



-oo 
oo 



-1^ / Pn{y)Pk{y)Mz) -y'o{y))w{y,t)dy 



oo 



— r- / Pn{y)Pk{y)wo{y)— — ^-]-dy 
hk J-oo y - 1 



oo 



-7^ / Pn{y)Pk{y){y'o{z) -yQ{y))w{y,t)dy 

l^k J-oo 

7 /"°° Pn{y)Pk{y) 



hk J-oo y - 1 



w{y,t)dy, (3.5) 



where we used notation \/o{z) = —\ogwo{z) . 

We note that analogous consideration of C„^„ = imphes 



pliyWoiyMy, t)dy = 7 / ^wiy, t)dy, (3.6) 

oo ^-oo y 

also the property C„ „_i = n implies the following Freud equation 

^ = T~I Pn{y)Pn-i{y)y'Q{y)w{y,t)dy- I ^""^^^^"'^^^^ w{y,t)dy. (3.7) 

Substitution of Cn,k into ()3.3jl and summation over k using the Christoffel-Darboux 
formula; 



Pj{x)Pj{y) _ Pn{x)Pn-l{y) - Pn{y)Pn-l{x) 



n-1 

y 

produces the differentiation formula; 

P'n{z) = -Bn{z)pn{z) + (3nAn{z)pn-l{z) , (3.^ 



where 



1 r v^,(^)-v'o(i/) 



An{z) := — / ''-^pi{y)w{y,t)dy + an{z,t) 

hn J-oo z-y 

7 p^iy) 

an{z,t) := — / ^ -w{y,t)dy 

hn J-oo {y-t){z-y) 

1 C"^ v' (z) - v' (v) 

Bn{z) := T / ^ '^Pn{y)Pn-l{y)w{y,t)dy + h^{z,t) 

hn-1 J-oo Z-y 

, / ,N 1 Pn{y)pn-i{y) I ,^ ^^ 

K[z,t) := / -w{y,t)dy 3.9 

hn-i J-oo {y-t){z-y) 
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Equation ()3.8|) is the "lowering" operator. 

A direct calculation produces two fundamental compatibility conditions valid for all z; 

Bn+l{z) + Bn{z) = {Z- an)An{z) - v'^iz) {Si) 
1 + {Z- an){Bn+l{z) - Bn{z)) = (3n+lAn+l{z) - /5„A„_i(z). (52) 

where we have used to arrive at (^i). Without going into details, we mention here that 
if the factor \x — t\'^ in the weight w{x,t) is replaced by 

AT 

then (5*1) and (5*2) still hold and the only changes are 

an{z,ti,..,t^) = J2h ( - M yM,-.tM)dy (3.10) 

j=i " J -00 [y '-jji^ y) 
hn{zM...M) = r !'''^^]\r"^^\ ^{yM....tu)dy. (3.11) 

^ K-i J -00 {y-tj)[z-y) 



(3.12) 

Using (5*1) and recurrence relations we have the "raising" operator, 

p'^_,iz) = (5„(z) + v'o(2;))p„_i(z) - Ar,.iiz)p^iz). (3.13) 

In the next section we take Wq{x) = exp(— x^), and make use of (Si) and (52) to produce 
a pair of non-linear difference equations satisfied by the recurrence coefficients for fixed t. 
These when combined with the t— evolution equations satisfied by the recurrence coefficients 
result in a particular Painleve IV. 

4 Derivation of the Painleve equation 

For w{x,t) = exp(— a;^)|x — t]'^, vq^x) = , we find, 

An{z) = 2 + an{z,t) (4.1) 

Bn{z) = bniz,t). (4.2) 

For z near 00, with fixed t, we obtain the following asymptotic expansions; 

/ ,^ 2an , 7 + 2ton , 7^ + 7«n + 2t2a„ 
an{z,t) ~ — H — \ — h .. (4.3) 

2(3n-n , t{2(3n-n) , ^ (3^ + t\2(3n - n) , 
On{Z,t) ~ \ \ h .., (4.4) 

z z^ z-^ 



where the coefficients are determined from orthogonahty, the recurrence relations, ()3.6p and 

(EZI). 

Substituting the asymptotic expansions into (^i) and (5*2), we find, by comparing the 
coefficients of ^ two difference equations satisfied by an and 

1 7 

(3n+i+Pn = n + - + - + an{t-an) (4.5) 

it - an) [j3n+l - Pn ~ ^ = Pn+l^n+l " /^nttn-l- (4.6) 

Remark 1. If 7 = 0, then a„ = 0, thus ()4.5|) and ()4.6|) become /?„+i+/5„ = n+1/2 and 
Pn+i—Pn = 1/2 respectively. The solution of these equations, subject to the initial condition 

= is Pn = which is the recurrence coefficients of the Hermite polynomials. 

Remark 2. If t = , then a„ = 0, then (4.5) becomes Pn+i + Pn = n + {1 + •y) /2. The 
unique solution subject to the initial condition /3o = is /5„ = n/2 + j{l — (— 1)")/4, which 
is the recurrence coefficient of what Szego called the generalized Hermite polynomials (see 
)27j . problem 25). These should not be confused with the generalized Hermite polynomials 
which arise in the rational solutions of Painleve IV (see next section). 

To study the t— evolution of the recurrence coefficients we begin by taking a derivative 
with respect to t of the squared norm hn of the n -th orthogonal polynomial, 

dthn 7 P PrXy) ( ,w o (A 7\ 

- -w{y,t)dy = -2an, (4.7) 



hn hn J ~oo t y 

where the last equality is obtained by using relation ()H.6j) and noting that y'^iy) = 2y. Since 
Pn = hn/hn-1, equatiou ()4.7|1 implies 

2ian-i-an). (4.8) 



Pn 
poo 

= / Pn{y)Pn-i{y)w{y,t)dy 



Differentiating relation 



with respect to t, we find, 

/oo 
Pn{y)Pn-i{y)wo{y)dt\y - tpdy, 
-00 

where function Pi{n) was defined in ()3.2|) . Using the Freud equation ()3.7p we now get 

^ / X 7 /"°° Pn{y)Pn-i{y) , ,^ , 
OtPlin) = / w[y,t)dy 



2 



h 



n-l 



00 



yPn{y)Pn~i{y)w{y, t)dy -n = 2pn-n. (4.9) 
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In view of relation ()3.ip . 

9ia„ = 2(/5„-/5„+i) + l. (4.10) 

The equations (j4.8p and (|4.1Up are the Toda evolution equations. 

We now show that D„ := D„exp(?T,t^) satisfies the Toda molecule equation (c.f. |231)- 
First note that ^ 

^aj = -Pi(n) and = -^"+^-^"-\ 
The equation (4.7) together with (4.9) implies 

n 



and hence 



a logD„ = 4- 



which is the Toda molecule equation. 

To proceed further, we parameterize as 



then relation (j4.5j) becomes 

^;^^±^ = (t -«„)«„. (4.12) 
Multiplying relation ()4.6|) by a„ and using the previous relation we get 

' n+1 n r) ^ 
1 = - Unan-lPn- 

Therefore 

where a does not depend on n. Taking into account the initial condition = — | we 
obtain the equation 

= 2 (n + r„ + -) + — . (4-13) 



In terms of the variables r„ , the Toda equations ()4.8|1 . ()4.1()j) become 



= an + 777 — -^^ — -—-dtTn, (4.14) 
2(n + r„ + I] 



and 

C^iftn = - ?^n+l (4.15) 



11 



respectively. Eliminating r^+i from equations ()4.15|) and (4.12) produces 

rn = oin{t - an) + ^<9tan- (4.16) 
To get the differential equation on we substitute expressions (I4.14j) . (j4.16p into (4.13): 

Oinit -an) + = (2n + 2a„(t - a„) + a'^ + 7)"^ + "n ( a„(t - a„) + ^ j + ^ 
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After simplification we obtain the following result. 

Theorem 1 The recurrent coefficients an{t) satisfy 

< = ^ + - ^tal + 2(^2 - 7 - 2n - l)a„ - ^ (4.17) 
which is a particular fourth Painleve equation. 

5 Explicit solutions for even multiplicity 

Painleve IV equation was first represented as a simple system of three first order equations 
(dressing chain) in [2H1- Such a symmetric form of PIV was used in [221 to obtain all the 
rational solutions of the equation in the remarkable determinant form (simultaneously with 
the independent work fE])- We use the notations from Noumi-Yamada to recall their 
results and then to use them. 

Firstly we bring equation ()4.17p to the canonical form by a simple change of variable. 
Let y = 2an and t = —t then ()4.17p takes the form 

y"it) = ^ + + ^ty' + 2(? - a)y + -, (5.1) 
2?/ 2 y 

where a = 2n + l + 7, b = —27^ . Then the symmetric form of PIV is a system of first order 
differential equations satisfied by /o = fo{x), fi = fi{x), f2 = f2{,x) , where 

/i(x) = -c|/(-cx), (5.2) 

with c = ^—3/2 . The system reads as follows: 

/o + /o(/i-/2) = bo (5.3) 

/{ + /i(/2-/o) = &i (5.4) 

/2 + /2(/o-/i) = h (5.5) 

where 

/0 + /1 + /2 = 3x (5.6) 

+ + = 3, (5.7) 
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and parameters of the PIV are suitably expressed in terms of bo, bi and 62- The PIV equation 
can also be written in the bilinear form on the level of r -functions. The solution of (j4.17|) 
may then be expressed in terms of r -functions tq{x) , Ti{x) , T2{x) as 

/i = -^log^ + x, (5.8) 

ax tq 

where the functions tq and T2 will be defined later. 

The generalised Hermite polynomials |22] are defined as 



Hm,ni^) = det (5.9) 

where 

i+2j=s •' 

They coincide with the specialization Sn^{x, |, 0, 0, . . .) of Schur polynomials corresponding 
to rectangular Young diagrams containing m rows of length n . 
Define also the set of functions 

/ \ f X^ ^ ~ ^ 2\ TT / \ 

Um,n[x) = ( 2 ^ j 

Then the triple 

i'^l^T2) ('^m.n; '^m+l,ru "^mjn+l) 

leads to a solution of PIV through formulas ()5.H) - ()5.8j) in the case 7 = m . 

Theorem 2 The recurrent coefficients an for the weight w{x) = e~^^ {x—tY^ with K G 
are given by 

Id H2K,n+i{t/c) 
an[t) = — — log 



2dt " H2K,nit/c) 



where Hm,n{x) are defined by L5. .9)) . and c = a/— 3/2 . 

Proof. For 7 = 2K with K eN the orthogonal polynomials with the weight w{x,t) can 
be expressed in terms of Hermite polynomials by the Christoffel formula ([2Z], pg. 30), since 
w{x,t) is the Hermite weight multiplied by a polynomials in x. It follows from the formula 
that the recurrence coefficients a„,/5„ are rational functions of t. Therefore an{t) is a 
rational solution of equation (j^.lj) in this case. The rational solution of the PIV equation is 
unique if it exists (see jT^I) and is expressed in terms of the generalized Hermite polynomials: 

(^n{t) = -y{-t) = -— /i - = -- — log- 



2c \cj 2dt H2K,n{t/c) 
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Remark. There is another way to see rationahty of an{t) . Indeed, equation ()2.13|) defines 
an even polynomials of degree 2Kn in t, hence, hn{t) = Dn+i{t)/Dn{t), is rational in t 
and (4.7) shows that an{t) is also rational in t. 

The above considerations allow us to obtain an expression for the Hankel determinant. 
We have seen that (jUTj), 

dthi 



IK 



Therefore 



So 



H2K,i+l{t C) 

log hi - log — — — = const. 

H2K,i+l{t/c) 

hi = ai— -— — 5.11 



for some constant a,- 



Proposition 2 (c.f. J^, The Hankel determinant for the weight w{x) = e ^^(x — t)" 

with K & is given by 

Dn = ^K,nH2K,n{tl C) 



where c = y— 3/2 , and 



A 



K-n 



||a, = (-l)^"7r2 



i=o 2^"+^^G(2K + 1) 

with the Barnes G— function defined by 

G{z + 1) = r{z)G{z), G(l) = l. 
Proof. It is clear from ()5.1H) and the product expression of Dn , 

Dn = h(jhi . . . hn-l 

where hi are the square of the norm of the monic orthogonal polynomials, that the 
constant A^^n in the proposition depends only on n,K , so all we need to do is to determine 
its value. 

Note that the coefficient of t^^" of Dn{t) is equal to the Hankel determinant associated 
with the Hermite weight. Therefore 

n— 1 .| 

Dn{t) = t^^'^-n^ n ^ + ^^'^^^ terms. (5.12) 

On the other hand the leading coefficient of H2K,n(t) is equal to 

G{2K + l)G{n + 1) 



G{2K + n+l) 
14 



(5.13) 



(see Combining ()5.12p and ()5.13p together we get the value of as stated. 

Remark. The Hankel determinant D„ as the average of characteristic polynomial (j2.13|) 
was first computed by Brezin and Hikami in ,5j as determinant of Hermite polynomials. 
The equivalence of the resulting formulas with the formulas for the r -functions H2K,n for 
PIV from [22] was used by Forrester and Witte in [12] (see also [IS])- We have now an 
explanation for this coincidence through showing that the diagonal recurrence coefficients 
a„(t) is a solution of PIV. We also note that this result can be obtained other way round 
using dlTj) and [T!T] . 
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